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A model to describe electronic correlations in energy bands is considered. The model is a gener- 
alization of the conventional Hubbard model that allows for the fact that the wavefunction for two 
electrons occupying the same Wannier orbital is different from the product of single electron wave- 
functions. We diagonalize the Hamiltonian exactly on a four-site cluster and study its properties as a 
function of band filling. The quasiparticle weight is found to decrease and the quasiparticle effective 
mass to increase as the electronic band filling increases, and spectral weight in one- and two-particle 
spectral functions is transfered from low to high frequencies as the band filling increases. Quasipar- 
ticles at the Fermi energy are found to be more 'dressed' when the Fermi level is in the upper half 
of the band (hole carriers) than when it is in the lower half of the band (electron carriers). The 
effective interaction between carriers is found to be strongly dependent on band filling becoming less 
repulsive as the band filling increases, and attractive near the top of the band in certain parameter 
ranges. The effective interaction is most attractive when the single hole carriers are most heavily 
dressed, and in the parameter regime where the effective interaction is attractive, hole carriers are 
found to 'undress', hence become more like electrons, when they pair. It is proposed that these are 
generic properties of electronic energy bands in solids that reflect a fundamental electron-hole asym- 
metry of condensed matter. The relation of these results to the understanding of superconductivity 
in solids is discussed. 
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I. INTRODUCTION 

The fact that the Coulomb repulsion between two elec- 
trons in a doubly occupied atomic orbital is smaller than 
the value predicted from the expectation value of the 
Coulomb interaction with single electron wavefunctions 
establishes that atomic orbitals are not infinitely rigid [Q: 
electrons will develop intra-atomic correlations to reduce 
their Coulomb repulsion. This well-known fact is not 
incorporated in the conventional single band Hubbard 
modelQ, which assumes that two electrons of opposite 
spin will doubly-occupy the same single electron atomic 
orbital. Dynamic Hubbard models || attempt to remedy 
this deficiency by either introducing an auxiliary boson 
degree of freedom to mimic the orbital relaxation that 
occurs on double atomic occupancy, or by allowing more 
than one atomic orbital. Here we consider a purely elec- 
tronic dynamic Hubbard model with two orbitals per site 
and no auxiliary boson degree of freedom. 

The model considered here has some superficial resem- 
blance to various multi-orbital Hubbard models that have 
been considered in the past such as the degenerate Hub- 
bard modelQ, the Falicov-Kimball model|| and the An- 
derson lattice model 0. However conceptually it is rather 
different. The goal is not to model the physics of elec- 
trons in degenerate atomic orbitals, nor of two different 
partially occupied bands in a solid, nor of local mag- 
netic moments interacting with conduction electrons, nor 
of mixed valence. Rather, we are interested in model- 
ing the physics of a single band, and based on general 



arguments Q argue that it is necessary, to understand 
single band physics, to include at least one other higher- 
lying orbital besides the one that is being filled as the 
electronic band is being filled. This second orbital be- 
comes increasingly important as the filling of the band 
increases beyond one half. 

When an electronic energy band is less than half-full, 
electrons can avoid each other by developing elaborate 
interatomic correlations. For this situation the ordinary 
single band Hubbard model may be adequate. However, 
when the band is more than half-filled, some Wannier 
orbitals are necessarily doubly occupied. The resulting 
large cost in Coulomb energy cannot be avoided by any 
interatomic correlations; instead, it can be reduced by 
the electrons developing intra-atomic correlations. The 
conventional single band Hubbard model cannot describe 
this, hence it is inadequate to describe electronic energy 
bands that are more than half-full: it is necessary to in- 
troduce another degree of freedom. Here this is achieved 
by having a second atomic orbital in the model, that be- 
comes increasingly occupied as the band filling increases 
beyond half filling. As we will see in this paper, the re- 
sulting electron-hole asymmetry in the band causes the 
quasiparticle weight of carriers in the upper half of the 
band (holes) to be smaller than that in the lower half 
of the band (electrons), and their effective mass to be 
larger. We believe that this physics is generally a part 
of the physics of electronic energy bands in solids. Its 
quantitative importance will depend on the particular 
solid under consideration. 

In section II we introduce the electronic dynamic Hub- 
bard model to be studied in this paper. Section III gives 
analytic results in limiting cases, and section IV presents 
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numerical results for effective interaction, quasiparticle 
weight and effective mass as a function of band filling 
and model parameters. Section V present results for op- 
tical absorption, and in Sect. VI we discuss the relation 
between the parameters in the Hamiltonian and atomic 
parameters in materials. We conclude in Sect. VII with 
a discussion of the applicability of these results to the 
description of real electrons in real solids. 

II. MODEL HAMILTONIAN 

The static (conventional) single band Hubbard model 
is defined by the Hamiltonian 

i,j,a i 

which properly emphasizes the fact that the intra-atomic 
electron-electron repulsion is the dominant source of elec- 
tronic correlation in solids. Here, cj CT creates an electron 
in a Wannier orbital centered at lattice site i, and the 
hopping amplitude Uj is the Fourier transform of the 
single electron band energy e^,: 

*« = ^£ ett( *" ,We *- ( 2 ) 

k 

The fundamental problem with the Hubbard Hamilto- 
nian Eq. (I) is that it implicitly assumes that the state 
of two electrons in a Wannier orbital is a single Slater 
determinant cJ^cJjO >, with |0 > the state of the empty 
orbital. This is incorrect, precisely because of the exis- 
tence of the electron-electron interaction described by the 
Hubbard U. We extend the Hamiltonian Eq. (1) to al- 
low for the fact that when two electrons occupy the same 
Wannier orbital, intra-orbital electronic correlations will 
develop due to the large local electron-electron repulsion. 
This effect can be described by having at least two Wan- 
nier orbitals per site, which the electrons will partially 
occupy to reduce their Coulomb repulsion. Let c| CT and 

c ia be the operators creating electrons into these Wannier 
orbitals. We consider the local Hamiltonian Q| 

Hi = Un^riii + U'n'^n'^ + en t + Vn^ (3a) 

a 

with n lry = c\ a c lrT , <r = c ia C 'ia> n * = 1*1 + n H> 
n\ = n'^ + n'^. In the extreme tight binding limit the 
Wannier orbitals become atomic orbitals, and we will cast 
the discussion in that language even though it should be 
more generally applicable. 

The primed orbital is higher in energy by e than the un- 
primed orbital, hence more extended in space, so that the 
Coulomb repulsion U' should be smaller than U. The or- 
bital energy e will generally satisfy the relation t < e < U, 
so that when the site is doubly occupied electrons will 



partially occupy this orbital to reduce their Coulomb re- 
pulsion and there is no justification for neglecting this 
second orbital as is done in the Hamiltonian Eq. (I). 
V describes the Coulomb repulsion between one electron 
in each orbital, and t' the intra-atomic hybridization be- 
tween these orbitals. The lattice Hamiltonian is then 

+ j> 

i 

with Hi given by Eq. (3a). We expect the various hop- 
ping matrix elements U, t^., t'^ to be similar in magnitude 
and hence will assume for the remainder of this paper 

t*j = tij = ^ij (4) 

and furthermore that the hopping connects only nearest 
neighbor sites, hence £y = t for i,j nearest neighbors, 
zero otherwise. 

It should be pointed out that we do not expect the 
qualitative physics of the model to depend on the par- 
ticular choice for the hoppings Eq. (4). For example, in 
Ref. § the same site Hamiltonian Eq. (3a) was consid- 
ered with hoppings t J • = ty = instead of Eq. (4) , and 
the same qualitative physics was obtained as in the model 
considered here. Certainly, quantitative differences will 
occur depending on these choices that should be studied 
in the future. In this connection it is also relevant to note 
that in a recent study of the periodic Anderson model, 
different choices of the hybridization (whether same site 
or neighboring sites) were found to not affect the under- 
lying physics [f§. 

III. ANALYTIC RESULTS 

A. Non-interacting band structure 

The single electron part of the Hamiltonian Eq. (3) is, 
after Fourier transforming 

« = E(iiHt) (-) 

with 

ffc = —2tcosk (6) 

for a one-dimensional geometry. Two bands result, with 
energy versus k relation 

E]: 2 =e k + e -±^{ e -)i + {e k -t>y (7) 
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The two bands are separated by an indirect gap 



A 



Ek=0 ~ E k=K 



(-y + (2t-t>y+ (8) 



which is always positive for e > 0. The bandwidth of 
the two bands is approximately given by Di = At + 2t', 
D 2 = 4i — 2t'. Figure 1 shows the band structure for 
various parameters. Note that the effective mass in the 
lower band is larger near the top of the band than near 
the bottom; this effect becomes much more pronounced 
in the presence of electron-electron interactions. 

The zero temperature single electron spectral function 
for the system with N electrons is 



where < Zk < 1 is the quasiparticle weight, and A' a is 
the incoherent part of the spectral function. We define 
the quasiparticle weight in our model at the Fermi energy 
by 

z(N) = \<0 N - 1 \c akFtT \0 N >\ 2 (15) 
with the quasiparticle operator defined by 



i i 

Cak F a — Uk F Ck F a T HfcpCj. a 



(16) 



A a { k ,w) 



I < n N+1 \cl k JO N > | 2 x 



with u\ + v\ = 1. We choose the particular linear combi- 
nation that maximizes the quasiparticle weight Eq. (15). 
The following simple relations are easily proven: 

#) = < O n ^\c 1<t \O n > | 2 + | < O n ^\cI\O n > | 5 

(17) 



5(w - (K 



JV+l 



N+l 




I < «w-i|c afc(7 |0jv > \ 2 S(u> + (E, 



N-l 



E 



N-l 



MAT 



(9) 
-i)) 



uL = E ' l<0 ^ l ^ >l (18a) 



z{N) 



Here, > | is the n-th excited state of the system 
with iV electrons (n = is the ground state) and E^ is 
the eigenvalue, and /ijy = E, 



N+l 



Eq , 



JV-1 



E, 



N 



Eq i . For a metal /ijv = (J-N-i = M an d we redefine the 
frequency u> — > lu + fi so that 

A a {k,u) = V| <n N+1 \cl ka \0 N > | 2 x (10) 



| < njv-ilcotolOjv > | 2 <5(^ + (E^ 1 - E*- 1 )) 

The index a labels the fcrmion operator of the unprimcd 
or primed orbital, or an appropriate linear combination 
theerof. The creation operator for an electron in band 1 
is given by 



Clka = U k C kcr + V k c' ko 



(11) 



with Uk, Vk obtained from diagonalization of Hk, Eq. 
(5b). Using this operator in Eq. (10) we have simply 



A 1 (k,u) = 8(w-(El-ri) 



(12) 



and the quasiparticle weight for a single electron in this 
band is the coefficient of the (5-function in Eq. (12) 



Zk = 1 



(13) 



B. Quasiparticle weight for interacting system 

In an interacting many-body system the spectral func- 
tion has the form 

A a (k, u) = z k S(w - (e k - //)) + A' a {k,u) (14) 



2 _ Y.i\<QN-iK a \o N >\ 2 , , 

z(N) (18b) 

Eq. (17) gives the total quasiparticle weight at the Fermi 
energy, and 1 — z gives the amount of spectral weight in 
the incoherent part of the spectral function. The quan- 
tities Uk F , Vk F indicate how much of the quasiparticle re- 
sides in the lower and upper orbital respectively. For a 
single band conventional Hubbard model the quasiparti- 
cle weight at the Fermi energy is given by the first term 
in Eq. (17). 



C. Optical conductivity and effective mass 

The current operator in our model is given by 

J = ^[(4+i<r + 4i CT )( c - + 4) - h.c] (19) 

i 

and we will compute the optical conductivity at zero tem- 
perature given by 

I < 01 Jim > I 2 „, , . , 

gl H=^E w w 6("-(E m -E )) (20) 

m Em ~ E ° 

The optical sum rule states that the integral of the optical 
conductivity is|| 

/>oo 

/ dwcriH = J < 0| -T\0 > (21) 
Jo 2 

with the kinetic energy given by 



T = ~* E( c t + 4)(<*+fe + (22) 

iad 
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In our two-orbital model, optical transitions include both 
'intraband' transitions as well as intcrband transitions to 
the second band. For parameters where a clear sepa- 
ration of energy scales occurs one can define an effective 
Hamiltonian describing the low energy part of the Hilbcrt 
space, and write a 'partial' conductivity sum rule 



dwo\ (u>) = Ai 



<0\-T eff \0> 



(23) 



where T e ff is the kinetic energy in the lower band, and 
the high frequency cutoff uj m excludes transitions to the 
second band. 



D. Strong coupling limit 

We consider the Hamiltonian Eq. (3) in the parameter 
regime 



U' + 2e < V + e < U 



U,U',V» e » t' 



(24a) 



(24b) 



These conditions ensure that a single electron at a site 
resides primarily in the lower orbital, while in the doubly 
occupied site two electrons of opposite spin reside pri- 
marily in the higher orbital. Some results in this limit 
were discussed in ref. [Q. The site eigenstates to lowest 
order in t' are 



>= | T> |0> +£|o > 1 1> 



(25a) 



IU >= |0 > | TI> +S'[\ T> |0 > +|0 >|> |0 >] (25b) 

where the first (second) ket in the product refers to the 
lower (higher) orbital, and 



S = t'/e 



5' 



V-U'-e 



(26a) 



(26b) 



The single site quasiparticle weight for a hole is, from Eq. 
(17) 

z h = \< INU > I 2 + I < Ilc'tlU >\ 2 = (S + 5'f = s 2 

(27a) 

and the single site quasiparticle weight for an electron is 



z e = |<0|c T |T>| 2 + |<0| C ' T |T>| 2 = l (27b) 

Hence the single site spectral function for a hole has other 
terms in addition to the quasiparticle term since z% < 
I, while the site spectral function for an electron is a 



single S— function. Consequently, in the solid the spectral 
function for a single electron in the ground state of the 
Hamiltonian Eq. (3) is a S— function, while the spectral 
function for a single hole in the lower band has incoherent 
contributions, as discussed in ref. ||. 

We can define quasiparticle operators for the lower 
band Ci„ as operators connecting the site ground states 



ct 



|0>= 



4io>=iu> 



(28a) 



(28b) 



and the relation between bare electron and quasiparticle 
operators isjnj 



[1 + (5" - l)rii,-a]ci a 



(29) 



Replacing the bare fermion operators in the kinetic en- 
ergy Eq. (3b) in terms of the quasiparticle operators 
yields an effective Hamiltonian for quasiparticles in the 
lower band, of the form 



H, 



eff 



-t^[l + (S-l)(n i ,_ CT +n J -_ CT )+ (30) 



(S - 1) h i ^ <J n j .^ a \{c\ cr Cj Cr + h.c.) 

This Hamiltonian gives rise to a 'correlated hopping' 
term which leads to pairing of holes. The enhancement of 
the hopping amplitude for a hole when the sites involved 
in the hopping process have an additional hole is 



At = tS(l - S) 



(31) 



and leads to superconductivity as discussed in detail 
elsewhere lfLIf. 



IV. NUMERICAL RESULTS 

We diagonalize the Hamiltonian Eq. (3) on an N = 4- 
site lattice with number of electrons N e ranging from to 
8, as appropriate to fill the lower band. The maximum 
number of states in the Hilbert space is 4900. To re- 
duce the computational complexity and because we will 
be interested in the case of large Coulomb interactions 
we discard all states where there are more than two elec- 
trons at any given site from the Hilbert space. We also 
discard states that have two electrons of the same spin at 
a site, because we will not be interested in spin-polarized 
states. This simplification should not change the quali- 
tative physics of interest here, and such states certainly 
become irrelevant in the limit of large interorbital repul- 
sion V; however it could have a quantitative effect for 
small V and should be investigated further. With these 
simplifications the maximum size of the Hilbert space is 
1024 states when there are 6 electrons in the four-site 
system. 
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Figure 2 shows schematically the band under consid- 
eration and the number of states in the Hilbcrt space for 
different positions of the Fermi level. Note that there are 
many more states in the Hilbert space for a given number 
of holes in the band than for the same number of elec- 
trons. This fact alone indicates that the system will be 
more incoherent when the band is more than half filled 
compared to when it is less than half filled. The figure 
shows one representative state in each sector: when there 
is one electron at the site, it resides primarily in the lower 
orbital, when there are two electrons at the site they are 
depicted as occupying the higher orbital which will be the 
dominant contribution in the parameter range of interest. 



A. Choice of boundary conditions 

For the small system under study boundary conditions 
are of course important. We have calculated properties 
of the system using periodic, antiperiodic and free ends 
boundary conditions (be). While periodic be may yield 
results closer to the thermodynamic limit in some cases, 
we believe free ends be are preferable for several reasons. 
Most importantly, the single particle eigenstates for the 
finite chain are non-degenerate with free ends be, while 
degenerate states occur for both periodic and antiperi- 
odic be. 

Consider the effective interaction between two elec- 
trons in the chain with N e electrons, defined as 

U eff (N e ) = E Q (N e ) + E (N e - 2) - 2E (N e - 1) (32) 

with Eo(N e ) the ground state energy with N e electrons. 
Figure 3 shows the effective interaction versus band occu- 
pation for various parameter values. For all the different 
boundary conditions the effective interaction becomes at- 
tractive near the top of the band for certain parameter 
values. This is a robust effect. On the other hand, it can 
be seen that the effective interaction can be negative for 
the half-filled band case for periodic be, and it is zero 
near the bottom of the band for antiperiodic be. These 
are spurious effects related to finite system degeneracies, 
that are expected to dissappear in the thermodynamic 
limit |l2]] . Instead, for free ends be the effective interac- 
tion is repulsive except near the top of the band. This be- 
havior should persist for larger systems. Quantitatively, 
we expect the magnitude of the attractive interaction to 
be somewhat smaller for the infinite system than for the 
small cluster [p"3[ . 

Another reason to use free ends boundary conditions is 
that the optical sum rule is satisfied in that case for the 
finite system |l4|], while if periodic boundary conditions 
are used the 'Drude weight' needs to be added by hand to 
the optical response. Furthermore, the expression for the 
quasiparticle weight Eq. (17) is not correct when there 
are degenerate states at the Fermi energy. Hence we will 
use free ends boundary conditions for the remainder of 
this paper. 



B. Quasiparticle bands 

We begin our study with a look at the energy level 
spectrum of the model. We choose sets of parameters 
that yield a clear separation of the spectrum of the lower 
band states. Figure 4 shows the energy levels as a func- 
tion of N e , the number of electrons in the cluster, for a 
non-interacting case (a) and for an interacting case (b). 
The dashed line separates the low-lying band 'intraband' 
states from the other states. For the noninteracting case, 
the spectrum in the lower band is nearly symmetric for 
electrons and holes, as one would expect. Instead, with 
electron-electron interactions a large electron-hole asym- 
metry exists. 

Consider first the states for a single electron (N e = 1) 
and for a single hole (N e — 7). There are 4 states in each 
case (= the number of sites in the cluster) and the dis- 
tance between the lowest and highest intraband state is 
the bandwidth. Clearly, for the interacting case the hole 
bandwidth is much smaller than the electron bandwidth. 
The effective mass of the quasiparticle is inversely propor- 
tional to the spacing between intraband states. Clearly, 
the quasihole is substantially heavier than the quasielec- 
tron in the presence of electron-electron interactions for 
these parameter values. 

As the number of electrons or of holes is increased, the 
number of intraband states increases. There are 6 intra- 
band states for the half-filled band (N e = 4), and 12 
states for all other occupations, because states with dou- 
ble site occupancy are pushed to much higher energies 
due to the large Coulomb repulsions. Note the asym- 
metry in the spectra in the interacting case for band 
filling less and more than one half: the quasiparticle 
band is substantially narrower when the carriers are holes 
(more than half-filled band) compared to the correspond- 
ing case when the carriers are electrons. In the effective 
strong coupling Hamiltonian one obtains a quasiparticle 
bandwidth that decreases monotonically as the electronic 
band filling increases: 

D(n e ) = D[l-^(l-S)r (33) 

with < n e = N e /N < 2 the band filling. 

C. Quasiparticle properties versus band filling 

We consider first a case where the second band is well 
separated in energy, with e = 10. Consider the evolution 
of the quasiparticle weight at the Fermi energy, given by 
Eq. (17), as the magnitude of the Coulomb interactions 
increase, shown in Fig. 5a. The quasiparticle weight is 1 
for the noninteracting case, and it remains 1 in the pres- 
ence of interactions at the bottom and at the top of the 
band, since the single electron and the single hole behave 
as free particles. When the band filling increases from 
empty or decreases from full the quasiparticle weight de- 
creases in the presence of interactions and is lowest at 



6 



the half filled band. Thus the spectral function will have 
largest incoherent contribution at and close to half-filling, 
as one expects in the conventional single band Hubbard 
model. 

The quasiparticlc weights in Fig. 5a appear to be 
electron-hole symmetric. However, a close look reveals 
that, except for the noninteracting case, the quasiparti- 
cle weight for holes is always slightly smaller than the 
corresponding one for electrons, i.e. 

z(n e ) > z{2 - n e ) (34) 

for n e < 1. This effect is due to the presence of the 
second band, and will exist always as long as e < oo in 
our model. 

Figure 5b shows the effective interaction defined by 
Eq. (32) versus band filling for these cases. The effective 
interaction becomes more repulsive as the bare repulsion 
parameters increase, and is approximately electron-hole 
symmetric in this case, as in the case of the conventional 
Hubbard model. 

We can estimate the effective mass or the effective hop- 
ping amplitude for the quasiparticle from the difference 
in energy between the ground state and the first excited 
state. Figure 5c shows the results for these cases. For 
the noninteracting case t e ff is approximately constant 
versus band filling, and as the repulsive interactions in- 
crease it decreases as the number of carriers in the band 
increases. Here the electron- hole asymmetry due to the 
fact that e < oo is more apparent, with holes always 
being heavier than electrons. 

Next we consider the effect of decreasing the inter- 
band energy separation e in the presence of Coulomb 
repulsion. Figure 6 shows results for e = 10, 5, 4 and 
2. As e decreases, the quasiparticle properties become 
increasingly electron-hole asymmetric. The quasiparti- 
cle weight, shown in Fig. 6a, is substantially smaller for 
holes than for electrons as e becomes small. Similarly the 
effective hopping amplitude shown in Fig. 6c decreases 
as the band filling increases, and holes are much heavier 
than electrons when e becomes small. The effective inter- 
action (Fig. 6b) remains repulsive for these parameters 
for all band fillings. 

For sufficiently small e and not too large U' however the 
effective interaction at the top of the band will become 
atractive. Figure 7 shows results for e = 2 and U' = 2. 
We also show for comparison the cases e = 2, U' = 5 and 
e = 5,U' = 2, where the effective interaction is always 
repulsive: both small e and small U' are required to yield 
an attractive interaction for holes. 

It is interesting to examine the quasiparticle weight 
and effective hopping amplitude for these cases. For 
the parameters where U e f / is attractive the quasiparticle 
weight for a single hole is smallest, and the quasiparticlc 
weight is larger for two holes (Fig 7a). This is in contrast 
to the cases of repulsive U e ff, where the quasiparticle 
weight for 2 holes is smaller than for 1 hole. Similarly 
the effective hopping (Figure 7c) for two holes is larger 
than for one hole when U e ff is attractive, and smaller 



when U e ff is repulsive. In other words, quasiparticles 
'undress', i.e. increase their quasiparticlc weight and de- 
crease their effective mass, when they pair. 

It is also interesting to examine the expectation value 
of the kinetic energy operator, Eq. (22). This is shown in 
Figure 8 for the three parameter sets under consideration. 
The kinetic energy is lowered both when electrons are 
added to the empty band and when holes are added to 
the full band. For the case of attractive U e ff where the 
single hole is highly dressed (small z) and the effective 
hopping is smallest, the kinetic energy is highest, as one 
would expect. As Fig. 8 shows, in that case only when a 
second hole is added to the full band the kinetic energy 
decreases below twice the value of the single hole kinetic 
energy. This indicates that pairing of holes is driven by 
lowering of kinetic energy in this model. 

Finally, Figure 9 shows the composition of the quasi- 
particles, as given by the quantities Uk F , Vk F in Eq. (18). 
As e decreases the quasiparticles occupy predominantly 
the higher orbital (large Vk f ) when the number of elec- 
trons in the band increases. In the case where the pair- 
ing interaction is attractive the quasiparticle weight for 
a hole in the system with 2 holes is also dominantly in 
the higher orbital, because of the large probability for 
two holes to be on the same site; in contrast, for small e 
and larger U' when the holes are not paired, Vk f is much 
smaller for the system with two holes because the holes 
occupy different sites. For the case of large e, the quasi- 
particle weight is dominantly in the lower orbital for all 
fillings. These results indicate that a necessary but not 
sufficient condition for pairing in this model is that pa- 
rameters are such that there is a large probability for 
electrons to occupy the higher orbital when the band is 
close to full. 



D. Hole pairing 

As seen in the previous section, the model can give rise 
to pairing for carriers near the top of the quasiparticle 
band with repulsive Coulomb interactions. In this section 
we examine in which regions of parameter space is the 
effective interaction between holes attractive, and shed 
some light on its origin. 

Figure 10 shows the dependence of U e ff for two holes 
at the top of the band on various Hamiltonian parame- 
ters. As seen in Fig. 10a, U e ff is attractive when the 
repulsive interaction in the higher orbital, U', is suffi- 
ciently small. This is because two electrons at a site 
will occupy dominantly the higher atomic orbital in the 
regime where pairing occurs. The energy difference be- 
tween the two atomic orbitals, e, plays an important role: 
both for large and for small e the attraction is suppressed, 
with e ~ 1.5 yielding the largest range of U' where U e /f 
is attractive. This value of e is found to be optimal for 
a wide range of the parameters U, V and t' , i.e. it is set 
by the value of the hopping amplitude t = 1. 

In contrast, attraction between holes is favored by a 
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large value of the lower orbital repulsion U, as seen in 
Figure 10b (we ignore the small region of very small U 
where U e f / is attractive which is presumably unphysical) . 
This is because the attraction requires a large change in 
the state of the remaining electron when a second elec- 
tron is removed from the site. If U becomes small, two 
electrons will occupy the smaller rather than the higher 
orbital and this effect is lost. Similarly, we find (not 
shown) that large values of the interorbital Coulomb re- 
pulsion V are favorable to pairing: the dependence of 
U e ff on V is similar to the dependence on U shown in 
Fig. 10b. 

As a function of the interorbital hybridization t' , pair- 
ing will occur when t' is not too large, as seen in Fig. 
10c. Again, the reason is presumably that the states of 
an electron in a singly and in a doubly occupied site need 
to be sufficiently different, which will not happen if the 
two orbitals are strongly mixed by t' . 

It is interesting to examine the change in kinetic energy 
when carriers pair. In Figure 11 we plot the difference 
between twice the kinetic energy of a hole in the filled 
band and that of two holes in the filled band: 

AT = 2 < T >lhole - < T >2holes (35) 

with the kinetic energy operator given by Eq. (22). It can 
be seen by comparison with Figure 10 that kinetic energy 
is always lowered (AT > 0) when carriers pair, i.e. in the 
regime where U e ff is attractive. The condition AT > 
is necessary but not sufficient to yield U e ff < 0. This is 
because pairing is associated with a decrease in kinetic 
energy and an increase in potential energy in this model. 

E. Comparison with conventional pairing 

In conventional models of superconductivity pairing 
arises from an effective electron-electron attraction in- 
duced by coupling to a boson degree of freedom that does 
not differentiate between electrons and holes. The result- 
ing effective interaction is electron-hole symmetric, and 
in such models pairing is driven by lowering of poten- 
tial rather than kinetic energy. We can describe such 
a scenario in our model by assuming negative values of 
the on-site interaction U, presumably resulting from in- 
tegrating out a boson. We also take a very large value 
of the interorbital spacing e so as to approach a single 
band Hubbard model, and compute the effective inter- 
action between 2 holes (which is the same as between 2 
electrons) from Eq. (32). Not surprisingly, U e ff for this 
conventional regime is attractive (repulsive) when U is 
attractive (repulsive) . 

It is interesting to compare the behavior of the model 
in the conventional regime (i.e. attractive Hubbard 
model regime) with that in the regime discussed in the 
previous subsection where holes pair, which we call 'dy- 
namic Hubbard' regime. We compute the effective inter- 
action U e ff versus U' in the dynamic Hubbard regime 
and versus U (positive and negative) in the conventional 



regimen, and in figure 12 we plot various properties as 
a function of the resulting effective interaction U e ff m 
both regimes. Figure 12a shows the ratio of quasiparti- 
cle weights in the system with two holes and one hole. In 
the conventional regime, zijz\ is less than 1 both when 
the effective interaction is attractive and when it is re- 
pulsive. In contrast, in the dynamic Hubbard regime the 
quasiparticle weight for the case of two holes becomes 
much larger than for a single hole when holes pair. Note 
however that Z2 can be bigger than z\ even when the 
effective interaction is still repulsive. In Fig. 12b we 
show the behavior of the effective hopping amplitude de- 
fined from the difference in energy between the ground 
state and first excited state in the model. In the conven- 
tional regime, holes becomes heavier than single particles 
{ti/ti < 1) when they pair, in contrast in the dynamic 
Hubbard regime paired holes are much lighter than sin- 
gle holes when pairing occurs. However once again pairs 
can be lighter than single holes even when the effective 
interaction is still repulsive. This is of course a finite size 
effect because in an infinite system for low hole concen- 
tration holes would be far from each other if not bound 
in a pair. Finally, Fig. 12c compares the change in ki- 
netic energy upon pairing. In the conventional regime the 
kinetic energy increases upon pairing, so that the differ- 
ence between twice the single hole kinetic energy and the 
pair kinetic energy is negative, while the kinetic energy 
decreases strongly upon pairing in the dynamic Hubbard 
regime. 

These results illustrate the qualitative difference in the 
physics of pairing in these two different regimes. In the 
regime of dynamic Hubbard physics, pairing is associ- 
ated with undressing |Tc|| , i.e. increase in quasiparticle 
weight, decrease in quasiparticle mass and lowering of 
kinetic energy. In the conventional regime, the physics 
of pairing is exactly opposite, pairing is associated with 
dressing, i.e. smaller quasiparticle weight, larger quasi- 
particle mass and increase in kinetic energy for the pair. 
In summary, upon pairing quasiparticles become more 
coherent and lighter in the dynamic Hubbard regime and 
more incoherent and heavier in the conventional regime. 

Note that the attractive Hubbard model can describe 
conventional superconductivity both in the short coher- 
ence length regime (large attractive U) and in the regime 
where the coherence length is thousands of lattice spac- 
ings (small attractive U). The qualitative contrast that 
we make here between the physics of that model and the 
physics of the dynamic Hubbard model applies to both 
regimes. 



V. OPTICAL CONDUCTIVITY 

We calculate the optical conductivity given by Eq. 
(20). The total optical spectral weight in the model is 
related to the expectation value of the kinetic energy op- 
erator as given by Eq. (21). We will divide the frequency 
range into a low frequency range with cutoff u> m = 2, 
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that defines the low frequency spectral weight Ai , and de- 
note the remaining spectral weight at higher frequencies 
by Ah- The low frequency spectral weight includes the 
'intra-band' spectral weight; it also includes some low fre- 
quency absorption that is not intra-band when the lower 
band is close to full and e is not too large. 

Figure 13 shows the dependence of the integrated op- 
tical spectral weights on band filling, for three sets of 
parameters. For Fig. 13(a), with e = 10, the absorption 
is approximately electron-hole symmetric; however even 
in this case with large e it can be seen that for holes the 
intra-band low frequency absorption is somewhat lower 
and the high frequency absorption is somewhat higher 
than for electrons. As e decreases (Fig. 13(b)) and even 
more so when U' also decreases (Fig. 13(c)), the intra- 
band absorption becomes much smaller for holes than 
for electrons. Note also that for the case of Fig. 13(b) 
where U e ff is still repulsive between holes the low fre- 
quency absorption for 2 holes is only slightly larger than 
for 1 hole; instead, as U' is decreased and U e ff becomes 
attractive (Fig. 13(c)) the intraband optical absorption 
for 2 holes becomes more than twice the intraband opti- 
cal absorption for 1 hole, because optical spectral weight 
is transfered from high to low frequencies when pairing 
occurs. 

Figure 14 compares the optical conductivity for the 
nearly empty and the nearly full band, for the cases with 
e = 10 and e = 2. For the large e case, the intraband con- 
ductivity (per particle) is only slightly smaller for holes 
than for electrons. Instead, for e — 2 there is a dra- 
matic difference in the optical conductivity for electrons 
and holes: for holes, the intraband conductivity is very 
small and most of the optical absorption occurs at higher 
frequencies. 

Next we consider the behavior of the optical conduc- 
tivity upon doping. For the case of large e, it is similar 
for electrons and for holes, as seen in Figure 15: the in- 
traband conductivity per carrier decreases slightly with 
doping, and some spectral weight is added at higher fre- 
quencies. For the case of small e, Figure 16, the behav- 
ior is similar for electrons but dramatically different for 
holes: in the latter case, there is a large increase in the 
low frequency spectral weight for the case of two holes, 
which is due to the undressing of holes when they pair. 
Furthermore there is an overall shift of the non-intraband 
spectral weight at higher frequencies to lower frequen- 
cies. Similar behavior is found in other realizations of 
dynamic Hubbard models jlj|. Such a transfer of opti- 
cal spectral weight from high to low frequencies has been 
seen in high T c cuprates upon hole doping and upon low- 
ering the temperature below the superconducting critical 
temperature Jl^, [l7| , |l§| , [Tg| , 
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Finally we show the behavior of optical absorption in 
the regime of conventional pairing. We choose a large 
value of the on-site attraction to illustrate the behav- 
ior clearly, however the qualitative behavior persists for 
smaller attractive interaction. The optical absorption is 
electron- hole symmetric as one would expect, and most 



of the optical spectral weight is at low frequencies (intra- 
band) in this case for all band fillings. Comparing the 
case of 1 hole and 2 holes (or 1 electron and 2 electrons) 
in Fig. 17 (b), it is seen that pairing is associated with 
a decrease in the low frequency optical spectral weight, 
i.e. quasiparticles become more dressed when they pair. 
This is in accordance with the behavior found for the 
quasiparticle weight and effective hopping in figure 12, 
and qualitatively different to the behavior in the dynamic 
Hubbard model regime. 



VI. RELATION WITH ATOMIC PHYSICS AND 
WITH REAL MATERIALS 

For any given atom one can relate the parameters in 
the site Hamiltonian 

Hi = U riii riii + U' n i-\ n 'ii + en i + V n i n 'i (36) 

a 

to atomic quantities by comparison of properties ob- 
tained from it and properties of the electronic states of 
the atom obtained from quantum chemical calculations. 
As the simplest example we discuss here qualitatively the 
relation between the Hamiltonian parameters and elec- 
trons in a hydrogenic ion with nuclear charge Z within 
the Hartree approximation. The difference in energy be- 
tween an electron in the Is and 2s atomic orbitals cor- 
responds to the energy difference between the two single 
particle eigenstates in Eq. (36), namely 



13.6 x Z 2 x 



(37) 



in eV units here and in what follows. For small t' we 
have approximately 



10.2^ 



(38) 



We will assume t' small in what follows so that the strong 
coupling analysis is applicable. The repulsion U in the 
lower orbital corresponds to the repulsion of 2 electrons 
in the Is orbital 



U = 17Z 



(39) 



In the Hartree approximation the single electron orbital 
with wavefunction tp oc e~ Zr expands upon double occu- 
pancy to wavefunction (p oc e~ Zr , with 



Z = Z 

16 



(40) 



We identify the Coulomb repulsion in the upper orbital, 
U', as the repulsion of two electrons in the Hartree ex- 
panded orbital, i.e. 



U' = 17Z = U- 5.31 



(41) 
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We can estimate the Coulomb repulsion between elec- 
trons in the two different orbitals, V, by calculating the 
Coulomb integral for one electron in the Is orbital and 
another in the expanded Hartree orbital. This yields (in 
eV) 



V = 27.2ZZ 



Z 2 + 3ZZ + Z 2 
{Z + Zf 



(42) 



Finally we can estimate the intra-atomic hopping t' from 
the overlap matrix element of the single electron wave 
function in the doubly and singly occupied sites. In the 
model, that is approximately 



s=<l\c\\n>=t'(\ + v -^ 7 

and in the Hartree atom it is 



S =< ip\ip >- 



(ZZ) 



3/2 



{{Z + Z)/2f 



(43) 



(44) 



It can be seen that in the atom U > V > U' for all Z. 
As Z decreases, all Coulomb repulsions decrease, as well 
as e and the overlap S. This is the regime favorable for 
pairing in this model. The Hartree calculation is of course 
very approximate, and in particular it overestimates the 
overlap matrix element 5. Nevertheless it illustrates the 
basic trend. For orbitals higher than the Is the energy 
levels become closer in energy and the effects discussed 
in this paper should become stronger. 

In summary, the atomic charge Z, with Z—2 the charge 
of the ion when the relevant band is full, is the key atomic 
parameter. For small Z the parameters in the Hamilto- 
nian studied in this paper move towards the regime of 
interest, namely small Coulomb repulsion U' , small in- 
terband separation e and small overlap matrix element 
S. In that regime electron-hole asymmetry in the band 
becomes dominant, holes become heavily dressed in the 
normal state and they strongly undress when they pair. 

For high T c cuprates the relevant band of interest is 
one formed by overlapping planar oxygen pn orbitals in 
the CuOi planes pH- Since in the undoped system (no 
holes) the ion is = , Z = in this case. For MgB2, the 
relevant band is formed by overlapping boron p xy orbitals 
in the B~ planes |||, and Z = 1. The fact that the 
planes are negatively charged in both cases (Z < 2) favors 
the physics discussed here, with the effects stronger for 
the cuprates due to the smaller Z. Even stronger hole 
dressing and higher T c 's would be expected in a structure 
with even smaller Z , for example if one managed to make 
a material with N = planes doped with some holes (Z = 

-!)■ 

The material LiBC has been recently proposed as a 
candidate for high temperature superconductivity when 
hole-doped, by analogy with MgB2, within electron- 
phonon theory |23||. Because the (BC)~ planes in that 
material would be less negatively charged than the (£?2) = 
in MgB^i i.e. effectively Z = 1.5 instead of Z = 1, we 
expect this not to be a modification of MgB2 conducive 



to higher T c 's within the physics discussed here. If such 
material was found to have a T c larger than MgB2, as 
predicted |2j|, it would directly contradict the assump- 
tions of this paper and prove the inapplicability of the 
concepts discussed here to real materials. 



VII. DISCUSSION 

Electrons in solids interact with each other with an in- 
teraction strength (e 2 = lAAeVA) that is of the same 
magnitude as the interaction strength of electrons with 
ions. It was recognized from the beginnings of solid state 
physics that Bloch's approach of prioritizing the electron- 
ion interaction over the electron-electron interaction was 
an ad-hoc assumption that could certainly not be rig- 
urously justified. Even though Landau's Fermi liquid 
theory with the concept of a quasiparticle provides an ex- 
planation for the fact that many properties of solids look 
amazingly 'independent-electron-like', the fundamental 
role of electron-electron interactions in solids is still not 
well understood. 

This paper is part of a continuing effort to understand 
the role of electronic correlation in energy bands. We 
argue that a key fact that has been ignored in previous 
treatments of the problem is the dependence of quasi- 
particle weight on band filling and the fundamental role 
of electron-hole asymmetry. In this paper we studied a 
'minimal model' that incorporates these key features. We 
believe that this physics is part of the physics of all elec- 
tronic energy bands: that quasiparticle weights at the 
Fermi level when the band filling n (0 < n < 2) is below 
and above the half- filled band (n = 1) are related by 



z(n) > z(2 - n) 



(45) 



with n < 1, i.e. that holes are always more dressed than 
electrons. How different the quasiparticle weights in the 
lower and upper halfs of the band arc determines how 
important the new physics originating in this effect is. 
This in turn depends on the ionic charge Z (Z — 2 =ionic 
charge when the band is full) with the strongest effects 
occurring for small Z . Because quite generally in an atom 
the intra-orbital Coulomb repulsion is linear in Z (e.g. 
Eq. (39)) and the energy level spacing is quadratic in 
Z (e.g. Eq. (38)) the effects discussed here will become 
unimportant for sufficiently large Z . 

We have called the models describing this physics 'dy- 
namic Hubbard models'. In these models, unlike the case 
of the conventional Hubbard model, the strength of the 
on-site repulsion U becomes a dynamical variable that 
can take more than one value depending on the state of 
the two electrons in the atom. Here this dynamics is 
incorporated by having two electronic orbitals per site; 
in other work we have described this dynamics with a 
single electronic orbital per site and an auxiliary boson 
degree of freedom jl3L 15, |24|, ^5|. While the model dis- 
cussed here is more realistic and closer to the physics of 
real atoms, the models with auxiliary boson degrees of 



10 



freedom are simpler to treat theoreticaly and thus may 
yield useful insight into the fundamental physics of this 
class of models. From the results in this paper and in 
previous work we believe that dynamic Hubbard models 
with only electronic degrees of freedom and those with 
auxiliary boson degrees of freedom share the same fun- 
damental physics. 

We have studied the two-orbital model by exact di- 
agonalization of a small cluster. It should be possi- 
ble to study larger clusters with more computing power 
and more sophisticated numerical techniques such as 
Lanczos diagonalization, density-matrix renormalization 
group and quantum Monte Carlo. We believe that the 
qualitative physics found here is likely to exist in larger 
systems. 

The calculations in this paper yield the properties of 
interacting electrons in a model Hamiltonian for the en- 
tire range of band fillings from empty to full, without un- 
controlled approximations. Before this work such studies 
had only been performed for simpler models such as the 
single band conventional Hubbard model, which as we 
have argued lacks some essential physics. The results 
found here should qualitatively apply to all electronic 
energy bands in solids. 

The results found here corroborate some of our ear- 
lier findings concerning the importance of electron-hole 
asymmetry fio|] and display clearly the interpolation be- 
tween the conventional understanding of electronic cor- 
relations in energy bands and the physics stressed in the 
theory of hole superconductivity. In the conventional un- 
derstanding electrons and holes are similar, quasiparti- 
cles are undressed when the band is almost empty and 
almost full, and the dressing and importance of electron- 
electron interactions increases as one approaches the half- 
filled band from either side. Instead, in the theory of 
hole superconductivity in its simplest interpretation the 
dressing of a quasiparticle increases monotonically as the 
band filling increases from the empty to the full band. 
As we have seen in this paper, the actual situation is al- 
ways in-between these two limiting descriptions, with the 
relationship Eq. (45) holding in all cases. 

The essential difference between conventional (static) 
and dynamic Hubbard models concerning 'intra-band' 
physics is that the state of a given electron is the same in 
the singly and doubly occupied atom in the static Hub- 
bard model, while it is different in the dynamic Hubbard 
model. Through this modification of the state the in- 
traband bare particles, which were strongly interacting 
with repulsion U, become weakly interacting quasipar- 
ticles with interaction U'. This occurs at the level of a 
single site, and is expressed by the relation between bare 
particle operators Ci a and quasiparticle operators Ci a 

Cia — [1 + (S - \)hi- a \ci a (46) 



The quasiparticle dynamics is described by the kinetic 
energy Eq. (30) and the local repulsion U', and their 
weight is further modified by the weak interactions in 
the quasiparticle band. The quasiparticle weight at the 
Fermi energy can be approximately written as 

z(n) = [l + (S-l)^z lb (n) (47) 

with < n < 2 the band filling, and the 'intraband' 
quasiparticle weight z^, defined by 

z t b = \<0N-i\c kF a\0N>\ 2 (48) 

calculated using the ground states of the Hamiltonian for 
the intraband quasiparticles (Eq. (30) plus weak on-site 
repulsion). In particular zu,(n — > 0) = z^in — ► 2) = 1 
and is smallest near the half-filled band, as in the conven- 
tional Hubbard model. The factor multiplying za in Eq. 
(47) isolates the main effect of electron-hole asymmetry. 
However even z^ will exhibit some additional electron- 
hole asymmetry (of the same sign) due to the dependence 
of the effective bandwidth on filling Eq. (33): the resid- 
ual intraband interactions will more strongly dress the 
quasiparticles in the upper half of the band where the 
effective bandwidth is smaller. 

For materials with large ionic charge Z , S will be close 
to 1, electrons and holes will be very similar and the dom- 
inant dressing will occur near the half-filled band. In- 
stead, for materials with small Z, S will be much smaller 
than 1, the physics of hole superconductivity will domi- 
nate, holes will be highly dressed near the full band and 
strongly undress as the local hole concentration increases. 
We propose that the physics of high temperature super- 
conductivity in solids is described by the latter regime. 
As the parameters become less extreme with increasing 
ionic charge Z the dressing of holes in the normal state 
becomes less extreme, the undressing effect of pairing be- 
comes less apparent, the coherence length of the Cooper 
pairs increases and one moves towards the regime of 'con- 
ventional' superconductivity pq|. 
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FIG. 2: Schematic depiction of states in 4-site cluster for dif- 
ferent band filling. Because of the large U in the lower orbital, 
two electrons at the same site will predominantly occupy the 
higher orbital. The number of states in the Hilbert space for 
each filling is also given. 



FIG. 3: Effective interaction Eq. (32) versus band filling 
for different boundary conditions, for four sets of parameters. 
t = 1 here and in the following figures. For all cases here 
U = 10 and V = 6. The lines through the data are a guide to 
the eye. Solid line: U' = 2, e = 2, t' = 0.2; dashed line: U' = 
4,e = 2,t' = 0.2; dash-dotted line: U' = 2, e = 1, t' = 0.2; 
dotted line: U' = 2, e = 2, t' = 0.5. 



FIG. 4: Energy eigenvalues for the different band fillings. Pa- 
rameters are: (a) e = 12, t' = 0.2, U = V = U' = 0; (b) 
e = 6, t' = 0.2, U = 20, V = 12, U' = 2. The dashed 
lines show the separation between the 'intra-band' states de- 
scribed approximately by the Hamiltonian Eq. (3) below it, 
and the rest of the states in the Hilbert space. For other less 
extreme parameters no clear separation is seen. Note that 
in the interacting case (b) the spectrum is strongly electron- 
hole asymmetric, with the lower band bandwidth for holes 
(N e > 4) much smaller than for electrons (N e < 4). 



FIG. 5: The parameters used are e = 10, t' = 0.2 and 
the following interactions: solid lines: U = V = U' = 0; 
dashed lines: U = 3, V = 2, U' = 1; dash-dotted lines: 
U = 6, V = 4, U' = 3; dotted lines: U = 10, V = 6, U' = 5. 
Plotted versus bandfilling N e are (a) quasiparticle weight at 
the Fermi energy, Eq. (17); (b) effective interaction Eq. (32); 
(c) effective hopping defined as the energy gap between the 
ground state and the first excited state. Note that even for 
this case of large e there is a small electron-hole asymme- 
try, with the quasiparticle weight and effective hopping being 
slightly smaller for holes than for electrons. 



FIG. 6: Same as figure 5 for t' = 0.2, U = 10, V = 6, U' = 5 
and e values given in the figures. As e decreases holes become 
more incoherent and heavier, i.e. smaller z and t e ff- 



FIG. 7: Same as figure 5 for t' = 0.2, U = 10, V = 6 and 
values of e and U' given in the figure. For small e and U' 
(solid lines) the effective interaction is attractive for holes; in 
that case, the quasiparticle weight at the Fermi energy for 
the band filled with 2 holes is larger than for the band with 1 
hole (a), and the effective hopping is larger for the band with 
2 holes than for the band with 1 hole. 



FIG. 8: Expectation value of the kinetic operator Eq. (22) 
for the same parameters as in Fig. 7. Note that for the case 
of attractive effective interaction (full line) the kinetic energy 
for two holes is lower than twice the kinetic energy for one 
hole, while in the cases of repulsive effective interaction it is 
higher. 



FIG. 9: Composition of the quasiparticles, from Eq. (18), for 
the parameters of Fig. 7. Uk F and Vk F give the amplitude of 
the quasiparticle at the Fermi energy in the lower and upper 
atomic orbital. As the band filling increases Uk F decreases 
and Vk F increases, with the changes being largest for the case 
of attractive effective interaction 

FIG. 10: Dependence of effective interaction for two holes in 
the filled band on Hamiltonian parameters. The values of e 
are given in the figures, (a) Versus U', with U = 10, V = 
6,t' = 0.2; (b) versus U, with U' = 2, V = 6, t' = 0.2; (c) 
versus t' , with U = 10, V = 6, U' = 2. 

FIG. 11: Dependence of kinetic energy difference between 
single holes and pair of holes, Eq. (35), on Hamiltonian pa- 
rameters, for the same cases as Fig. 9. 

FIG. 12: Comparison of behavior of model in conventional 
(electron-hole symmetric) regime and regime of dynamic Hub- 
bard physics. Parameters used for dynamic Hubbard regime 
are U = 10, V = 6, e = 2, t' = 0.2 and U' ranging from to 6; 
for conventional regime, e = 100, t' — 0.2, V = 6, U' = 5 and 
U ranging from -1.5 to 8. In both cases f/ e // is calculated 
from Eq. (32) and the results are plotted versus U e ff. (a) 
Ratio of quasiparticle weights for two holes and for one hole; 
(b) ratio of hopping amplitudes for two holes and one hole; (c) 
difference in kinetic energy of single holes and paired holes, 
Eq. (35). 

FIG. 13: Integrals of optical conductivity. Ai denotes the 
low frequency integral Eq. (23), with cutoff uj m = 2; Ah 
denotes the high frequency optical spectral weight for uo > u) m , 
and A — Ai + Ah the total optical spectral weight Eq. (21). 
Parameters are U = 10, V = 6, t' = 0.2 and values of U' and 
e given in the figures (same parameters as Figs. 7-9). 

FIG. 14: Comparison of optical conductivity for one electron 
and one hole for t' — 0.2 and values of e given in the figure. 
Here and in the following figures of optical conductivity, the <5- 
functions in Eq. (20) are broadened to Lorentzians with width 
r = 0.25. The lowest frequency S— function at the 'Drude 
precursor' frequency is shifted to u> — and represented 
by a Drude form (semi-Lorentzian) . Note that for large e 
the conductivities for electrons and holes are very similar, 
for small e the conductivity for holes is very small and in 
particular the intra-band conductivity is much smaller than 
for the case of large e. 

FIG. 15: Comparison of optical conductivity for 1 carrier and 
2 carriers in the band. The optical conductivity for one carrier 
is multiplied by a factor 2 to make it comparable to the optical 
conductivity for two carriers, (a) electrons, (b) holes. Here 
and in the next figure U = 10, V = 6, t' = 0.2. Here, e = 
10, U' = 5. Note that the optical conductivity is similar for 
1 and 2 carriers (normalized to number of carriers) both for 
electrons and for holes, with the low frequency absorption (per 
carrier) being slightly smaller when the number of carriers is 
larger. 
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FIG. 16: Same as Fig. 15 for parameters e = 2, U' = 2. For 
the case of electrons the results are similar to Fig. 15, for the 
case of holes they are very different: for the case of 2 holes 
there is a large increase in low frequency absorption and an 
overall shift in optical absorption to lower frequencies. 



FIG. 17: Optical absorption for parameters in the 'conven- 
tional' regime, given in the caption of Fig. 12, with strong 
on-site attractive interaction U = —8, giving rise to effective 
attraction U e ff = —5.6. (a) Integrated optical absorption, 
same cutoff as in Fig. 13. Note that the low frequency ab- 
sorption does not increase as the number of carriers increases 
from 1 to 2. (b) Optical conductivity for 1 and 2 holes (es- 
sentially the same as for 1 and 2 electrons for these param- 
eters). Note that upon pairing the low frequency absorption 
decreases strongly as paired carriers are more highly dressed. 
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